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ABSTRACT 

We use the semi-analytic model GallCS to predict the Tully-Fisher relation in the B, I 
and K bands, and its evolution with redshift, up to z ~ 1. We refined the determination 
of the disk galaxies rotation velocity, with a dynamical recipe for the rotation curve, 
rather than a simple conversion from the total mass to maximum velocity. The new 
recipe takes into account the disk shape factor, and the angular momentum transfer 
occurring during secular evolution leading to the formation of bulges. This produces 
model rotation velocities that are lower by ~ 30 km/s in case of Milky Way-like objects, 
and ^ 20 — 30 km/s for the majority of the spirals, amounting to an average effect of 
~ 20 — 25%. We implemented stellar population models with a complete treatment of 
the thermally pulsing asymptotic giant branch, which leads to a revision of the mass- 
to-light ratio in the near-IR. Due to this effect, K band luminosities increase by ^ 0.5 
at redshift z = and by > 1 at z = 3, while in the I band at the same redshifts the 
increase amounts to ^ 0.3 and ^ 0.5 mags. With these two new recipes in place, the 
comparison between the predicted Tully-Fisher relation with a series of datasets in the 
optical and near- infrared, at redshifts between and 1, is used as a diagnostics of the 
assembly and evolution of spiral galaxies in the model. At redshifts 0.4 < z < 1.2 the 
match between the new model and data is remarkably good, expecially for later-type 
spirals (Sb/Sc). At z = the new model shows a net improvement in comparison with 
its original version of 2003, and in accord with recent observations in the K band, the 
model Tully-Fisher also shows a morphological differentiation. However, in all bands 
the z = model Tully-Fisher is too bright. We argue that this behaviour is caused by 
inadequate star formation histories in the model galaxies at low redshifts. The star- 
formation rate declines too slowly, due to continuous gas infall that is not efficiently 
suppressed. An analysis of the model disk scale lengths, at odds with observations, 
hints to some missing physics in the modeling of disk formation inside dark matter 
halos. 

Key words: galaxies: formation galaxies: evolution galaxies: kinematics and dynam- 
ics galaxies: structure galaxies: fundamental parameters galaxies: photometry 



1 INTRODUCTION 

The theory of galaxy formation has advanced remarkably in 
the past decade, thanks to the employment of sophisticated 
methods of investigation such as numerical simulations and 
semi-analytic models (see for instance Balland et al. 2003, 
Baugh et al. 2005, Bower et al. 2006, Cattaneo et al. 2008, 
Cole et al. 2000, De Lucia et al. 2004, Firmani & Avila-Reese 
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1999, Hatton et al. 2003, Kauffmann et al. 1993, Menci et 
al. 2006, Monaco et al. 2007, Somerville et al. 2008, Trujillo- 
Gomez et al. 2010, van den Bosch 2000, van den Bosch 2002). 
With these tools, the problem of the assembly of structures 
in the Universe can be addressed, and the vast number of 
non-linear physical processes that accompany the formation 
and evolution of galaxies from primordial perturbations can 
be modelled. 

Interesting for galaxy formation are the observations of 
scaling relations, that trace a pattern in the assembly his- 
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tory of the vast variety of objects that we observe in the 
Universe. The origin and nature of such relations is at the 
core of the mechanism of galaxy assembly (Courteau et al. 
2007). One of the most firmly established scaling relations is 
the TuUy-Fishcr (hereafter TF; TuUy & Fisher, 1977), that 
correlates the luminosity and the rotation velocity of disk 
galaxies. The tightness and universality of the relation make 
it a valid tool to measure extragalactic luminosity distances, 
since the rotational velocity does not depend on distance or 
cosmology. But a most interesting aspect of the relation is 
also the insight it provides on the mechanism of disk assem- 
bly (Tonini et al. 2006a). 

In the cold dark matter (CDM) scenario for disk galaxy 
formation (Fall & Efstathiou, 1980), the disk forms out of 
gas that cools inside a dark matter halo. Initially in equi- 
librium with the dark matter, the gas collapses into a ro- 
tationally supported disk under angular momentum conser- 
vation, while the halo evolves adiabatically (Blumcnthal et 
al., 1986). As a consequence, the final equilibrium structure 
and dynamics of the disk are strongly correlated with those 
of the dark matter halo, thus providing an insight into the 
physical halo parameters. In particular, the disk rotation 
velocity Vb, the disk size parameterized through the scale- 
length Rr>, and the disk mass Mo are determined by the 
halo virial properties, Vvii oc i?vir oc M^/f (Button et al., 
2007). The slope of the relation is a function of the radius 
at which the velocity is measured, mirroring the dark mat- 
ter mass distribution as a function of galactocentric distance 
(Yegorova & Salucci, 2007). 

Interestingly, although both dynamics and baryonic 
physics enter the TF, the relation is remarkably tight. In 
particular, galaxies follow the baryonic TF, which links the 
dynamical properties of disk (stars + gas) and dark matter 
halo, with a very small scatter (see for instance McGaugh 
et al. 2009). In the luminosity vs velocity TF relation, the 
measured luminosity is given by the stellar omission, and is 
therefore a proxy for the stellar component. The star forma- 
tion rate (SFR) and the feedback processes, that regulate 
the balance between stellar and gas mass, follow a scaling 
with the dynamical quantities, and provide at the same time 
a source of scatter. 

In galaxy formation models, the supernova feedback ef- 
ficiency is one of the main parameters driving the slope and 
zero-point of the TF relation. In fact, being the circular 
velocity of the disk linked to the halo virial velocity, and 
therefore to its escape velocity, it is also closely related to 
the ability of the galax;y to blow a wind. The balance be- 
tween feedback and star formation sets the amplitude of the 
TF relation, and its variation across the mass range affects 
the TF slope (Croton et al. 2006, dc Lucia ot al. 2004). 

Semi-analytic models so far do not have a good grasp of 
the TF relation. In fact, although both Croton et al. (2006) 
and de Lucia et al. (2004) reproduce the local I-band TF 
relation (Giovannelli et al. 1997), this match should not be 
considered a success because they do not model rotation 
curves in a physically justified way. De Lucia et al., assumes 
the total mass profiles are isothermal (i.e., constant circular 
velocity), while Croton et al. assume the observed rotation 
velocity is equal to the iiiaxinmm circular velocity of the 
dark matter halo (in the absense of galax;y formation) . Both 
of these assumptions are expected to under-predict the true 
rotation velocities. 



On the other hand, semi-analytic models which actu- 
ally model the galaxy rotation curves do not match the TF 
relation. Among these, GallCS (Hatton et al. 2003) includes 
baryons. Cole et al. (2000), Benson et al. (2003) and Ben- 
son & Bower (2010) include baryons, NFW halos (Navarro, 
Frenk & White (1997) and halo contraction. 

In addition, it has been pointed out that semi-analytic 
models cannot simultaneously account for the mass and 
angular momentum distribution of disk galaxies in a gas- 
dynamical context (Courteau et al., 2007, van den Bosch ct 
al. 2002b, Governato et al. 2007). In particular for the mod- 
eling of the TF relation, sources of systematic errors other 
than the feedback implementation include the structure of 
dark matter halos inferred from simulations (the halo con- 
centration and density profile for instance), and the model- 
ing of the dynamics of cooling and baryonic collapse inside 
halos (see for instance Tonini et al. 2006a, 2006b, Piontek 
& Steinmetz 2009). 

Of interest in the present paper, the predicted TF rela- 
tion is also affected by the modeling of the galaxy rotational 
velocity, and the stellar population models implemented. 
In this work we improve both these recipes in the GallCS 
model, and rc-addrcss the comparison with the observed TF 
relation. We consider a broad spectral range from the optical 
to the near-IR, and in particular we perform the compari- 
son in the K band. We also compare the model and observed 
evolution of the TF up to redshifts « ~ 1, both in the optical 
and in the near-IR. 

In current semi-analytic models, the rotation velocity 
of disk galaxies is a crude estimate based on the virial ve- 
locity of the host dark matter halos (for instance de Lucia 
et al. 2004, Croton ct al. 2006), or the sum of a halo -I- stel- 
lar component in spherical symmetry (Hatton et al. 2003). 
These approximations do not take into account the actual 
mass profile of the disk, the effects of angular momentum 
exchange in the galaxy, and the dark matter halo expansion 
due to dynamical friction, and introduce systematic errors 
both in the zero-point and in the slope of the predicted TF 
(see Button et al. 2007). In support to this point, Monaco 
et al. (2007) also showed that the observed baryonic TF 
relation can be matched by hierarchical models if the dark 
matter halo concentration is lower than predicted by ACBM 
simulations (to c = 4 for halo masses around 10^^ Ms). 

In the present work, we implement for the first time 
in the semi-analytic model GallCS (Hatton et al. 2003) a 
recipe for the galaxy rotation curve that takes into account 
the shape of the disk and the angular momentum exchange 
between disk and bulge and between the galaxy and the 
dark matter halo. We do so by modifying the disk and halo 
parameters following the recipe proposed by Button et al. 
(2007). 

The TF relation predicted by galaxy formation models 
is also determined by the conversion between stellar mass 
and light, which is performed through the use of stellar pop- 
ulation models. In the original version of GallCS (Hatton et 
al. 2003) the observed TF relation was not reproduced. In 
two recent papers (Tonini et al. 2009, 2010) we showed how 
the predictions of semi-analytic models are greatly affected 
by the choice of the input stellar populations. In particu- 
lar we showed how the use of the Maraston (2005; hereafter 
M05) models allow GallCS to reproduce the colours and 
near-IR luminosities of 2: ~ 2 galaxies, in the mass and age 
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ranges accessible to hierarchical clustering. In this paper we 
address the impact of the M05 models on the predicted TF 
relation in the optical and near-IR, and its evolution with 
redshift. 

This Paper is organized as follows. In Section 2 we in- 
troduce the galaxy formation model, and we describe the 
theoretical recipe to derive the galaxy rotation curves from 
the dynamical properties of the galaxies and their host dark 
matter halos, the stellar population models employed, and 
the method for the selection of model spirals, including the 
definition of a proxy for the morphological differentiation of 
these objects. In Section 3 we describe the data samples used 
to make the comparison between model and observations. In 
Section 4 we describe the results for the z — Tully-Fisher, 
and its dependence on morphology. In Section 5 we address 
the evolution of the TF with redshift. Finally, in Section 6 
we discuss our results and conclude. 



2 THE MODEL 

We produce the model galax;ies through the hybrid semi- 
analytic model GallCS (Hatton ct al. 2003), and we defer 
the reader to its original paper for details on the dark matter 
N-body simulation and the implementation of the baryonic 
physics. In brief, the model builds up the galaxies hierar- 
chically, and evolves the metallicity consistently with the 
cooling and star formation history (with the new implemen- 
tation of the chemical evolution by Pipino et al., 2009). Feed- 
back recipes for supernovae-driven winds and AGN activity 
are implemented in the code (the latest with the improved 
version of Cattaneo et al. 2008). Merger-driven morphology 
evolution and satellite stripping and disruption are taken 
into account. 

In order to compare the model output with the data, 
we produce mock galaxy catalogues both in rest-frame and 
observer-frame, adopting each time the broadband filter 
sets of the data available for comparison. The broadband 
magnitudes thus obtained are further scattered with gaus- 
sian errors comparable to the observational errors of our 
data samples (on average cr = 0.1 mag). Dust extinction 
is treated in the model, however the observed TF relations 
that we use for the comparison with data are corrected for 
internal extinction by the various authors. For this reason, 
we produce unreddened galaxies for the comparison with 
rest-frame TF relations. When comparing our model galax- 
ies with observed-frame data (not corrected for reddening), 
we implement dust extinction with a Calzetti extinction 
curve and a colour-excess E{B — V) proportional to the 
star formation rate for each single galaxy, parameterized as 
E{B -V) = 1/3 • (iogiSFR) - 2) + 1/3 (see Tonini et al. 
2010 for a discussion about this choice). Magnitudes are pre- 
sented in the AB system. We calculate the galaxy rotational 
velocity at the standard value of r = 2.2Ro (where Ru is 
the exponential disk scale-length), consistently with the ob- 
servations used in our comparisons. 

The semi-analytic model is rendered for the parameter 
set Ho = 71.9 km s"^ Mpc"\ Oa = 0.742 and Qm = 0.258; 
these parameters arc used to compute obscrvcd-framc mag- 
nitudes. The model absolute magnitudes are not directly 
dependent on the cosmology (which obviously enters in the 
physics of the semi-analytic model, but not in the compu- 



tation of the rest-frame light emission from galaxies); when 
comparing with multiple sets of data, analyzed by different 
authors with various cosmologies, we re-scale all units to the 
model h = i?o/100. 



2.1 Galaxy rotational velocity 

We calculate the rotational velocity of each galaxy by build- 
ing its rotation curve and evaluating it at r — 2.2Ro, where 
Ro is the disk scale-length. Galaxies in the models are sys- 
tems composed of a dark matter halo and stars distributed 

in a disk and a bulge; neglecting the gas component (which 
is not dominant in mass, at least at z < 1) the rotation curve 
is 



Vrltir) = Vi{r) + V-d'isk(r) + Vb'uige(r-) , 



(1) 



with r being the galactocentric distance. For the spheroidal 
components, i.e. the halo and the bulge, the contribution to 
the velocity is simply V^{r) = GM{r)/r. Following the 
original GallCS recipe, the halo density profile is defined as 
a truncated isothermal sphere. In the simulation halos are 
tidally truncated by the surrounding density field, and the 
inner region of halos is virialized. For the purpose of building 
the rotation curves, we consider halos as spherically symmet- 
ric isothermal spheres out to the virial radius. In addition, 
halos have a finite density core to avoid the singularity. The 
density profile can be parameterized as follows: 



p(r) = 



Po 



l/a 



:«{/3-7) 



(2) 



with (a,/3,7) = (1,2,0). Fits to the radial mass distrubu- 
tion of halos in dark matter-only N-body simulations pro- 
duce the NFW density profile (NFW 1997), which belongs 
to the same family of functions, with parameters (1, 3, 1). At 
intermediate radii NFW halos and isothermal spheres have 
the same density profile, while at the centre the NFW is 
shallower and at the outskirts it is steeper. However, in the 
range of galactocentric distances of interest in the present 
work, the difference between the two profiles is negligible 
(as is illustrated later on in this Section). In general, the 
choice of the halo functional form is a source of uncertainty, 
and it should be noted that there is an ongoing debate in 
the literature regarding the shape of dark matter halos in 
the presence of baryons, expecially in the center of galaxies 
where the dark matter component reacts to baryonic mass 
assembly with either expansion or contraction (see for in- 
stance Tonini et al. 2006b, Duffy et al. 2010). 

The halo is almost entirely pressure-supported, with a 
mild spin component, parameterized by A = J|£|^/^/GM^''^ 
in terms of the halo total mass M, angular momentum J 
and binding energy E. In simulations the distribution of the 
halo spin parameter appears to be log-normal, with both 
median and scatter independent of mass and redshift to a 
good approximation (see e.g. Bullock et al. 2001, Vitvitska 
et al. 2002, Hetznecker & Burkert 2006, Bett et al. 2007, 
Maccio et al. 2007). The median spin parameter for GallCS 
relaxed halos is A ~ 0.04. For each galaxy in GallCS, all halo 
parameters (virial mass and radius, and velocity dispersions) 
are directly taken from the N-body simulation, and naturally 
cause some scatter in the halo profile. 
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The bulge stellar mass distribution follows a Hernquist 
(1990) profile, characterized by parameters (1,4,1), and is 
assumed to be entirely pressure-supported (see Hatton et al. 
2003 for details on bulge formation and geometry). 

The disk structure and rotation is determined by the 
halo spin parameter A following Mo et al. (1998). The disk 
forms from a baryonic gaseous component that initially has 
the same specific angular momentum as the halo, a mass mo 
in units of the halo mass, and a total angular momentum jd 
in units of the halo angular momentum. The disk is assumed 
to have an exponential surface density profile, with a scale- 
length defined as: 




where -Rh is the halo virial radius (see Mo et al. 1998 for 
details). In galaxy formation models the assumption of spe- 
cific angular momentum conservation during disk formation 
is widely used, which leads to j-o/m-a = 1. This assump- 
tion, while useful to semplificate the problem and avoid ad- 
ditional degrees of freedom, is often unrealistic. In fact, an- 
gular momentum transfer from the infalling baryons to the 
halo pushes the ratio lower than unity, and both the halo and 
the disk structure are be affected (see Tonini et al. 2006b; 
also, Piontek & Steinmetz 2009). On the other hand, the 
ratio can also become greater than unity in case of removal 
of low angular momentum material (e.g.. Mailer & Dekel 
2002, Button & van den Bosch 2009). A discussion on the 
disk sizes of the model galaxies is to be found in the last 
Section. 

In the Hatton et al. (2003) version of GallCS, the 
disk component of the velocity was treated as V^{r) = 
GM{r) /r (for future reference, this will be called the 
'spheroidal' velocity model). However, neglecting the shape 
factor in the disk gravitational potential leads to an un- 
derestimation of the maximum rotational velocity by about 
15%, under angular momentum conservation (Binney & 
Tremaine, 2008). For the disk component, in order to take 
into account the shape factor we assume that 

VL^ir) = 1^ (-^) ' [loKo - hK^l/2R^ , (4) 

where Mo is the total disk mass, and /, K are modified 
Bessel functions of first and second kind (Tonini et al. 2006, 
Salucci et al. 2007). This recipe effectively increases the ro- 
tational velocity of disk-dominated systems, and has there- 
fore an effect on the TF relation expecially at the high-mass 
end. For future reference, we call Vdisky the total rotational 
velocity when V^j^^y. is modeled as in Eq. Q. 

Angular momentum conservation is not a realistic sce- 
nario also during bulge formation, as pointed out by Button 
et al. (2007). When bulges form through disk instabilities, 
part of the material in the disk loses its angular momen- 
tum, leading to the formation of a pressure-supported com- 
ponent at the center. The lost angular momentum is trans- 
ferred to the rest of the disk and the dark matter halo, as 
seen in numerical simulations (Hohl 1971, Bebattista et al. 
2006). This affects the value of Ro through jn/rno after 
bulge formation. Following Button et al. (2007), after defin- 
ing = A/buige/A^disk, the ucw disk scale-length is 

Ru = -^M /h[1 + (1 - ,fx)e] Ai?H ; (5) 
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Figure 1. Rotation curve of a simulated galaxy at z = 0, with 
disk mass M^isk = 3.9 lO"*^" Mq and bulge mass M^uigc = 
3.2 10^0 Mq. The 'new', 'disky' and 'spheroidal' models for the 
rotation curve are compared (from upper to lower panel). The 
vertical dotted green line shows the radius r = 2.2i?o before in- 
cluding the angular momentum transfer from bulge to disk. In the 
upper panel, the vertical solid green line shows the new location of 
r = 2.2iJD after angular momentum transfer. In all panels, solid 
thick black lines show the total rotation curve, solid blue lines 
show the disk component of the velocity, solid red lines show the 
bulge component. The solid thin black lines show the halo com- 
ponent as a truncated isothermal sphere, while the dashed thin 
black lines show the halo contribution as an NFW profile. The 
dashed thick black lines show the total rotation curve if the halo 
is an NFW. 



the quantity /x expresses the ratio between the specific an- 
gular momentum lost by the disk to the halo during bulge 
formation, to the total specific angular momentum of the 
baryonic material that went to form disk and bulge. We 
re-calculate the disk scale-lengths of the galaxies in the sim- 
ulation with this new recipe, obtaining 

i?D_ = iiDMMW (l + (1 - -M^M^S) ' (6) 

and we use the fiducial value /x = 0.25 indicated by Button 
et al. (2007). By inserting this new value of Ro into Eq. Q, 
we obtain what we will refer to as the 'new' model for the 
total rotational velocity, V^icw, from Eq. ([l|. 

Fig. ((T)) shows the comparison between the three models 
for the rotation curve, 'new', 'disky' and 'spheroidal' {panels 
from upper to lower, respectively), for a simulated galaxy at 
z — 0, which was chosen for the conspicuous bulge (Afbuigo ~ 
3.2 10^0 Mq and Mdisk = 3.9 10^° Mq). In all panels, the 
vertical dotted green line shows the radius r = 2.2Ro before 
the evaluation of the angular momentum transfer from bulge 
to disk. In the upper panel, the new value of r = 2.2Ro 
following the angular momentum transfer ('new' model) is 
shown as a vertical solid green line. In all panels, solid thick 
black lines show the total rotation curve, solid blue lines 
show the disk component of the velocity, solid red lines show 
the bulge component. 
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Note that, for a pure disk, the peak of the rotation 
curve is always approximately at r = 2.2Ro (making this 
particular radius the favourite in observations). From the 
'spheroidal' to the 'disky' model, the disk component of the 
velocity features a more distinguished peak, with the result 
that the maximum of the total rotation curve is closer to 
r = 2.2Ro, and the bulge is less prominent in determining its 
value. The sharpening of the disk peak causes the measured 
total velocity to increase. 

However, the presence of other mass components, such 
as the bulge and the halo, tend to broaden the total veloc- 
ity peak and push it to smaller radii, where these compo- 
nents are more concentrated. But the concentration of the 
bulge mass at the expense of the disk must be accompa- 
nied by angular momentum transfer from the bulge to the 
disk itself, and according to Eq. ((6]), the higher the ratio 
Mbuigo/Mdisk, the more the disk scale-lenght consequently 
increases, leading to the 'new' model and a re-distribution 
of the disk mass. The disk expands radially and, since it 
conserves its mass, it becomes less dense. Its velocity peak 
broadens and it is pushed outwards, thus stretching the to- 
tal rotation curve. Due to the decreased density, the disk 
potential well is lower, and this, together with the rapidly 
declining bulge contribution, lowers the total velocity mea- 
sured at the new r = 2.2Ro. Since this effect if brought 
about by the presence of the bulge, it is going to be ex- 
pecially important at the high-mass end of the model TF 
relation, where bulge formation is more frequent (this is a 
feature of the hierarchical model which will be discussed 
later on). 

In Fig. ([1} we also show the halo component of the total 
rotation curve. The thin solid black line is the velocity profile 
of the truncated isothermal sphere, while the thin dashed 
black line is the velocity profile of the NFW halo of same 
virial parameters. The total rotation curve obtained with the 
NFW {thick solid black lines) is virtually undistinguishable 
from the one obtained with the truncated isothermal sphere 
{solid thick black lines), expecially around the region r = 
2.2Ru and beyond. 

The new recipe for the rotation velocity allows for a 
sharper differentiation of the model rotation curves depend- 
ing on the galaxy morphology, a fact that has been long since 
observed in nature (Salucci et al. 2007). In particular, in the 
upper panel of Fig. ((TJ the bump in the rotation curve due 
to the presence of the bulge component is clearly visible. 

The velocity recipe can be further optimized if one takes 
into account the dark matter halo evolution triggered by the 
formation of the galaxy. As proposed in Tonini et al. (2006b) 
and discussed in Button et al. (2007), angular momentum 
exchange between the collapsing baryonic component and 
the dark matter causes halo expansion, which affects the 
halo structural shape and ultimately the dark matter con- 
tribution to the rotation curve. However, such a recipe would 
require more fundamental modifications to the prescriptions 
for the dark matter halos in the semi-analytic code, and is 
beyond the scope of this paper. 

Fig. ([21 illustrates the difference in the predicted total 
rotational velocity (at 2 = 0) at the radius r = 2.2Ru be- 
tween the 'new' velocity model and the 'spheroidal' model 
{magenta dots), and between the 'new' model and the 'disky' 
model {cyan dots), as a function of log{W = 2Vncw), which 
is used for the comparison with observations later on. 



spilcroidal new 
disky new 
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Figure 2. The difference in the predicted 2 = galaxy rotational 
velocity at r = 2.2/Jd with different assumptions for the disk ve- 
locity, as a function of the total 'new' rotation velocity. Magenta 
points: difference between the original Hatton et al. (2003) model 
and the 'new' model. Cyan points): difference between the 'disky' 
model and the 'new' model presented here. For each of the mod- 
els, we plot the total velocity difference AV against the quantity 
log{W = 2Vrot), which is used for the comparison with observa- 
tions later on. 

The combined effect of increasing the disk scale-lenght 
and decreasing the bulge contribution to the rotation curve, 
is that the 'new' model velocity is always lower than the 
'spheroidal' model velocity, with a difference up to ~ 30 
km/s at the high-mass end. On the other hand, the 'new' 
model velocity is also always lower than the 'disky' model 
velocity, by up to ~ 50 km/s at the high-mass end, where 
bulges contribute significantly to the stellar mass. The dif- 
ference AV correlates in fact strongly with the total mass, 
and tends to zero for galaxies with small bulges. 

With the new recipe in place, model rotation veloci- 
ties are lower by ~ 20 — 30 km/s or less, compared to the 
'spheroidal' model, for spirals less massive than the Milky 
Way, and by ~ 40 — 50 km/s for Milky Way- like objects. 

In what follows, we adopt the 'new' model to compare 
the model predictions with the data. 

2.2 Stellar population models 

The photometry of the model galaxies is produced with the 
use of the M05 stellar population models, which include a 
detailed treatment of the Thermally Pulsing Asymptotic Gi- 
ant Branch phase of stellar evolution. The implementation 
of the M05 stellar populations in the semi-analytic model 
significantly improves the predictions for the colours and 
near-IR luminosities of galaxies, expecially at high redshift. 
The comparison of the model performance when equipped 
with the M05 and with its original input stellar popula- 
tions (PEGASE, Fioc & Rocca-Volmerange, 1997) shows 
that the inclusion of the TP-AGB produces near-IR lumi- 
nosities more than 1 magnitude brighter for a given stellar 
mass, and V — K colours more than 1 magnitude redder, 
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Figure 3. Difference between the predicted rest-frame TF in the 
M05 and PEGASE runs, as a function of redshift. The green dots 
show the I band and the red dots show the K band. 



matching the observations of high-redshift galaxies (Tonini 
et ah 2009, 2010, Henriques et aL in prep.). 

The TP-AGB hght emission starts at wavelengths A > 
5000 A, peaks in the rest-frame K band, and affects the 
nearby bands from the red part of the optical spectrum 
down to the infrared. This affects both the zero-point and 
the slope of the predicted TF relation, as a function of red- 
shift. Fig. ((3)l shows the difference in magnitude between the 
TF relation predicted with the M05 and the PEGASE runs 
of the semi-analytic model, in the rest-frame I band {green 
dots) and K band {red dots), in four redshift bins from 
to 3. On the a;— axis we plot W = 2Vto%- The difference be- 
tween the M05 and PEGASE TF in the I band tends to 
mildly increase with increasing redshift. The median magni- 
tude difference in the redshift bins centered in 2 = 0, 1, 2, 3 is 
respectively A(TF) ~ 0.15,0.22,0.25,0.2 mags. At z = 2, 3 
there is a mild dependence with galaxy mass, so that the dif- 
ference can go up to A(r_F) ~ 0.35—0.4 for massive galaxies. 
In the K band, we find a significant increase of the differ- 
ence A(r_F) with redshift. In the redshift bins centered in 
2 = 0, 1, 2, 3 we find the median A(TF) ~ 0.2, 0.65, 0.8, 0.7 
mags with maximum A(rF) > 1 mags for 2^1. There is 
no significant trend with galaxy mass. 

From now on in this work, the M05 models will be used 
for all the comparison with data. 



2.3 Selection of model galaxies 

The determination of the morphology of the model galaxies 
is based on the relative contribution of bulge and disk to 
some galactic property, and can be carried about in various 
ways. For the purpose of a more comprehensive study, we 
adopted two different criteria for the selection of the model 
spiral galaxies. 

Hatton et al. (2003) defined as spiral galaxies the ob- 
jects for which the ratio between the bulge and disk lu- 
minosity in the B band did not exceed a given threshold. 



Simien & de Vaucouleurs (1986) found that this ratio cor- 
relates well with the Hubble type (see for instance Graham 
& Worley 2008 for an updated study). Following this cri- 
terium, if we define the parameter 7 = e2;j3(— Lbuigo/idisk), 
the objects with 7 > 0.507 are defined as spirals (see also 
Baugh et al. 1996). This classification has the advantage 
to be observationally-oriented, but it depends on the spec- 
trophotometric model in use (including dust reddening), and 
yields no dynamical information. 

From the semi-analytic model point of view, we have 
the information about the mass and size of disk and bulge 
before we compute the stellar emission, so we can also clas- 
sify our galaxies directly in terms of the ratio between the 
bulge and disk stellar mass, to mimic the morphological dif- 
ferentiation. We choose to split our galaxies in 5 subsam- 
ples, defined by Mbuigo/Mdisk in 5 intervals: galaxies with 
Afbuigo/Mdiak < 0.25 loosely correspond to Sc/Irregulars, 
galaxies with 0.25 < Mbuigc/Afdisk < 0.5 loosely correspond 
to Sb spirals, galaxies with 0.5 < Afbuigo/Afdisk < 0.75 
loosely correspond to Sa spirals, and galaxies with 0.75 < 
Afbuigo/Afdisk < 1 are similar to 'SO' types. Strictly speak- 
ing, observed SO galaxies are not identified depending on the 
bulge/disk mass ratio, but are rather selected to be passively 
evolving disks, based on colours. Our morphology selection 
does not depend on photometry, but nonetheless it naturally 
preserves the correspondence between bulges and old stellar 
populations, since bulges in the model do not form new stars 
(except in merger events). 

We find that our somewhat arbitrary selection criterium 
qualitatively matches the observed properties of spirals (K. 
Masters, private communication). It is also consistent with 
other prescriptions used in the literature: we re-selected our 
Sb/Sc spirals following both Croton et al. (2006) and de Lu- 
cia et al. (2004), and found that the resulting subsample co- 
incides with the one selected in terms of Mbuigc/Mdisk- The 
totality of the spirals in the range < Afbuigc/Miisk < 1 is 
almost coincident with the sample selected with the Hatton 
et al. (2003) criterium (with the exception of some massive 
bulge-dominated objects that Hatton selected as spirals) . As 
Afbuigc/A'/disk increases, the average age of the stellar popu- 
lations tends naturally to increase. 

The hierarchical nature of the semi-analytic model 
causes a large number of spiral galaxies to be subject to 
some kind of perturbations (like gas stripping and strangu- 
lation) when they live in dense regions, with the consequent 
quenching of star formation. In any galaxy identified as a 
satellite, i.e. infalling into the halo of a central galaxy (which 
on the contrary sits in the center of its halo), cooling of fresh 
gas is shut down, unless received via a merger. A compari- 
son of the model satellites with observed spirals is therefore 
improper, and as will be shown in Section 4, largely unsuc- 
cessful. For this reason, we use only model central galaxies 
to build the TF. 



3 DATA SELECTION 

The data at our disposal are either in rest-frame or observed- 
frame. When data are in rest-frame, the absolute magni- 
tudes, evolutionary and fc-corrections, and corrections for 
internal extinction were performed by the various authors. 
We compare these sets of data with unreddened model galax- 
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ies. When data are in observed-frame and uncorrected for 
dust extinction, we compare them with model galaxies to 
which we apply our reddening recipe. 



3.1 Local data 

For the z — TuUy-Fisher relation in the I and K bands, 
wc use the data by Masters ct al. (2006, 2008). Their I-band 
sample, called Spiral Field I++ is a compilation of objects 
taken from the Spiral Field I band and the Spiral Cluster I 
band catalogues analyzed by Giovanelli et al. (1994, 1995, 
1997a,b) and Hayncs ct al. (1999a,b), the SC2 sample (Dale 
1998, Dale et al. 1999), data from Vogt (1995) and Catinella 
(2005) and the HI archive from Springbob et al. (2005). The 
SFI++ contains 5000 galaxies, and Masters et al. (2006) an- 
alyze a subsample of data consisting of spirals in the vicinity 
of nearby clusters. The rotational widths arc derived from HI 
global profiles when possible (about 60% of the sample), or 
optical rotation curves. See Masters et al. (2006) for details 
about the velocity derivation. These objects are specifically 
targeted for I-band observations, and although the sample is 
not complete (in magnitude or volume), it includes all types 
of spirals. However, as a result of various bias corrections, 
the sample is designed to be dominated by late-type spirals 
(see Masters et al. 2006). The derived TF relation is consis- 
tent with the one published by Giovannclli ct al. (1997). 

For the K band, the data come from the 2MASS Tully- 
Fisher Survey (2MTF), cross-matched with the 2MASS Ex- 
tended Source Catalog (XSC) including all the galaxies with 
rotational widths coming either from the Cornell HI digi- 
tal archive (Springbob et al. 2005) or the Cornell database 
of Optical Rotation Curves (ORCs) described by Catinella 
(2005). The cross-match was performed for objects in the 
vicinity of the clusters in the SFI-I— |- sample. About 65% of 
the galaxies in the 2MTF sample belong also to the SFI++ 
sample. Details about the corrections applied to the ro- 
tational widths, to account for inclination, turbulent mo- 
tions, cosmological broadening and instrumental effects can 
be found in Masters et al. (2008), along with a detailed list 
of all the bias corrections. 

For both the I and K band samples described above. 
Masters et al. (2006,2008) find a morphological dependence 
of the slope and zero-point of the Tully-Fisher relation. 

Further data for the I band, together with the B band 
at z = 0, come from Verheijen (2001). In this case the sam- 
ple belongs to the Ursa Major cluster (TuUy et al. 1996). 
Only galaxies with a lower inclination limit of 45° are con- 
sidered, which leaves a sample of 49 galaxies described as 
the complete sample. Corrections for internal extinction are 
treated following TuUy ct al. (1998). In the complete sam- 
ple, the rotation velocity is estimated from the width of HI 
lines, corrected for instrumental resolution, internal turbu- 
lent motions and inclination. Further subsamplcs arc defined 
according to the observed properties of the galaxies; one of 
these subsamples is called the Rotation Curve sample (RC), 
containing galaxies for which synthesis imaging data is avail- 
able, and the shape of the HI rotation curve is measured. In 
this case, velocity is also estimated at the maximum of the 
rotation curve (Knax) or in the flat part of the curve (Vflat) 
depending on the curve shape; these estimates are in good 
agreement with the Hl-width method (see Verheijen 2001 



for details). We also use the relation derived in the B band 
by Pierce & Tully (1992) for comparison. 

3.2 High-redshift data 

The evolution of the B-band TF relation up to ~ 1.3 is 
compared with data from Fernandez-Lorenzo et al. (2009), 

in a sample belonging to the DEEP2 survey (Davis et al., 
2003). The photometric catalogue (Data Release 1, Coil et 
al. 2004) features magnitudes corrected for dust extinction 
(after Tully et al. 1998), and is ft-corrected following Blan- 
ton & Roweis (2007). The rotation velocity has been es- 
timated from optical line widths of integrated spectra, ob- 
tained through line-fitting techniques. The selection of spiral 
galaxies was performed through visual inspection of HST 
images of the AEGIS survey, where it overlaps with the 
DEEP2. Inclination angles of the selected sample are all be- 
tween 25° and 80°. The sample was purged of E/SO galaxies 
and interactive pairs, and includes spirals and irregulars. For 
comparison, we also take the published TF of Bamford et 
al. (2006) for field galaxies at < z >= 0.33. 

The K-band TF up to z ~ 1.2 is compared with data 
from Fernandez-Lorenzo ct al. (2010). The sample belongs 
to the Groth Strip Survey, with spectroscopy provided by 
the DEEP2 survey (Data Release 3). The photometric data 
arc part the AEGIS survey (Davis et al. 2003, 2007); the K- 
band data were collected with Palomar WIRC (Bundy et al. 
2006), and with the Infrared Array Camera on the Spitzer 
Space Telescope (Barmby et al. 2008). Data are corrected 
for internal extinction, and galaxy morphology is selected 
by eye on Hubble Space Telescope images. As in Fernandez- 
Lorenzo et al. (2009), the rotational velocity is determined 
through optical line widths. 

Another sample of optical and near-IR data comes from 
Bohm & Ziegler (2007). It contains 73 galaxies with pho- 
tometry in 7 bands from the UV to the near-IR. This ob- 
servational sample was constructed utilizing the multi-band 
imaging survey of the FORS Deep Field (FDF, Heidt et al. 
2003), which comprises deep U, B, g, R, I, J and K pho- 
tometry performed with the Very Large Telescope and the 
New Technology Telescope. To extract rotation curves from 
the 2D spectra, the A-positions (linecenters) of a given emis- 
sion line were measured by fitting Gaussians stepwise along 
its spatial profile. Due to the small apparent sizes of the 
galaxies, the observed rotation curves are heavily blurred. 
At redshifts oi z ~ 0.5, the apparent scale lengths are of the 
same order as the slit width and the seeing disk. To correct 
for these blurring efi^ects, synthetic rotation velocity fields 
were generated assuming an intrinsic shape with a linear 
rise of the rotation speed at small radii and a turnover into 
a regime of constant rotation velocity Vmax at large radii 
(for details on this prescription as well as tests of various 
rotation curve shapes, see Bohm et al. 2004). Intrinsic dust 
absorption was corrected following Tully & Fouque (1985). 



4 Z=0 TULLY-FISHER RELATION 

In this Section wc present the 2 = model TF compared 
to the observational relations in the B, I and K photometric 
bands. This choice of bands is functional to test the overall 
performance of the model in relation to the star formation 
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rates and the balance between stellar populations of difTer- 
ent ages in our galaxies, as well as the contribution of the 
TP-AGB light in the near-IR. We show the morphological 
differentiation of the K-band TF relation, which is predicted 
in the model and present in the data. Moreover, we analyze 
the relative contributions of our dynamical recipe and the 
stellar population models in determining the K-band TF re- 
lation. We compare the performance of the model for central 
galaxies and satellites. We also discuss a possible origin for 
SO galaxies. 

4.1 Tully-Fisher relation in B, I and K bands 

Figure Q shows the predicted z = Q Tully-Fisher rela- 
tion in the B, I and K bands [from left to right); scat- 
ter plots are shown in the upper panels, while the median 
relations are shown in the lower panels. Disk galaxies in 
the model are selected according to the Mbuige/Mdisk cri- 
terium: yellow for A/buigc /A/disk < 0.25, magenta for 0.25 < 
Mbuige/Mdisk < 0.5, blue for 0.5 < Mbuige/Afdisk < 0.75 and 
red for 0.75 < Afbuigc/Afdisk < 1. The whole sample of spiral 
galaxies selected following Hatton et al. (2003) in terms of 
the relative luminosity of bulge and disk is represented by 
the cyan dots. The model velocities are taken at the fiducial 
value of 2.2 Ru- In all bands, we compare the model pre- 
dictions with the observational relation by Verheijen (2001), 
which is represented by two lines, indicating a different se- 
lection criterium for the observed gala:xies. The green solid 
line) represents the complete sample with a measured HI 
global profile, while the green dashed line) represent the RC 
subsample. In the B band, we also plot the observational re- 
lation from Pierce & TuUy (1992; thick solid black line). In 
the I and K bands, we compare the model TF with the obser- 
vational relation by Masters et al. (2006, 2008, black solid 
line for the relation, and dot-dashed lines for its scatter). 
The I-band TF relation by Masters et al. (2006) is coinci- 
dent with the relation found by Giovannelli et al. (1997). 

For the B band, the plot shows that the model is suc- 
cessful in reproducing the observed TF slope, expecially for 
the later-type disk galaxies, and lies instead on the bright 
side in terms of the zero point. The offset with the Pierce & 
Tully relation is substantial, while it amounts to ~ 0.2 — 0.4 
mags with respect to the Verheijen relation. Notice the dis- 
crepancy between different observational relations as well; 
in particular, the relation by Pierce & Tully (1992) is on 
the faint side of the others, and in fact led to a measured 
value of the Hubble parameter Ho ~ 85 km/s/Mpc, much 
larger than more recent results. For this reason, we consider 
the discrepancy between the model TF and this particular 
observational relation less significant. 

In the I band, the model is able to reproduce the slope 
of the I-band TF relation very well, and is offset in the zero- 
point by ~ 0.2 — 0.4 mags. The model scatter is also com- 
parable to the observed one. The performance of the model 
TF is worse here than the results of Croton et al. (2006) 
and de Lucia et al. (2004), who compare their models with 
the observed relation of Giovannelli et al. (1997). Both the 
Croton et al. (2006) and de Lucia et al. (2004) models differ 
from ours in the stellar population models adopted; by us- 
ing BC03 instead of M05, the TF is shifted towards fainter 
magnitudes. Moreover, their oversimplification of the rota- 
tion velocity recipe shifts the TF towards slower values. 



In the K band, the agreement of the model predictions 
with the data is better than in the other bands. Consistently, 
the K band is less affected by uncertainties such as dust red- 
dening. Although the model relation is accomodated inside 
the data scatter, we nonetheless notice a discrepancy in the 
slope, which is more evident than in other bands. We can 
also identify a cut in velocity, around log{W) ~ 2.3, below 
which the model galaxies are too bright. 

Table lU lists the slopes and zero-points at logW = 2.4 
of the model TF relations shown in Fig. Q, obtained with 
power-law fits of the kind M — 51ogh = a -I- blogW, sepa- 
rated by morphological types. Interestingly, in all the bands 
considered the model shows a net differentiation between 
morphological types, with the more disk-dominated galaxies 
following a steeper relation with brighter zero-points in all 
bands, while bulge-dominated objects show a flatter slope. 
The morphology dependece of the K-band TF will be inves- 
tigated in detail later on. 

Notice the turn in the model TF slope at the high-mass 
end, for earlier types; this is the combined effect of AGN 
feedback and bulge formation. On the one hand, in the more 
massive objects the onset of AGN activity quenches star 
formation and thus reduces the galaxy luminosity per unit 
mass, moving the objects below the relation. On the other 
hand, while at the low-mass end the mechanism of bulge 
formation is preferentially secular evolution, which does not 
introduce a large scatter in the galaxy stellar populations, at 
the high-mass end bulges form primarily through mergers, 
which contribute to increase the scatter and raise the mass- 
to-light ratio. In addition, the larger bulge masses boost the 
bulge contribution to the rotation curve and increases the 
rotation velocity. These factors combined push the model 
galaxies below the relation at the high-mass end, and in- 
crease the scatter in the relation. This may also be the cause 
why the velocity range covered by GallCS spirals is not as 
wide as the data. The model in fact does not seem to produce 
late-type galaxies that are massive enough (contrary to both 
the De Lucia et al. 2004 and Croton et al. 2006 models). If 
GallCS overpredicts the abundance of large bulges at the 
high-mass end, then the number of massive spirals evolving 
into early-types increases, and these objects are lost from the 
spiral selection. This is mirrored in the progressively larger 
scatter and flattening of the slope for earlier-type model spi- 
rals. 

Overall, although the model produces a reasonable TF 
slope, its performance is unsatisfactory in the B band, it is 
marginally acceptable in the I band, and it is better, but not 
perfect, in the K band. The interplay between dynamics and 
baryonic physics that shapes the TF relation is affected by 
many factors. We argue that we have good control over the 
dynamical model for the galaxy rotation curves and the stel- 
lar population models, both of which are physically sensible. 
Other degrees of freedom are introduced by the supernovae 
feedback recipe, the chemical evolution model, the star for- 
mation history of the model galaxies, dust extiction, and 
the modeling of the data to obtain rest-frame quantities (al- 
though such adjustments as the fc-correction are small for 
local samples). The increasingly good match between model 
and data from the optical to the near-infrared suggests two 
possible causes for the offset. 

The first possible explanation is the treatment of dust 
reddening. The observed relations used for the comparison 
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Figure 4. The 2: = Tully-Fisher relation in the B, I and K bands {from left to right). The upper panels show the scatter plots, while 
the lower panels show the median relations. In all panels, model spiral galaxies are selected according to Mbuigc/A^digj^ (see text). (Left 
panel): the B-band model TF {left panel) is compared with the observed relation by Verheijen (2001, green solid/dashed lines, see text) 
and Pierce & TuUy (1992; solid black line). {Central panel): the I-band model TF is compared with the observed one from Masters 
et al. 2006 {black solid line for the relation, and dot-dashed lines for its scatter) and Verheijen 2001 {green solid/dashed lines). {Right 
panel): the K-band model TF relation is compared with the observed relation by Masters et al. (2008; solid black line for the relation, 
and dot-dashed lines for its scatter) and Verheijen 2001 {green solid/dashed lines). 



Table 1.^ = model TF relations, power-law fits 



[A-^bulge/A'-^disk] 


[0,0.25] 


[0.25,0.5] 


[0.5,0.75] 


[0.75,1.] 


all spirals 


B band 


slope AAf/AlogW 


-7.48 


-7.12 


-6.52 


-5.20 


-4.95 


zero point at logW = 2.4 


-20.76 


-20.48 


-20.15 


-19.91 


-20.18 


I band 


slope AA//AlogW 


-7.51 


-6.66 


-6.43 


-5.33 


-5.37 


zero point at logW = 2.4 


-21.50 


-21.33 


-21.12 


-20.94 


-21.11 


K band 


slope AA//AlogW 


-7.97 


-7.02 


-6.74 


-5.72 


-5.71 


zero point at logW = 2.4 


-21.68 


-21.47 


-21.26 


-21.07 


-21.26 



have been corrected for internal extinction. Verheijen (2001) 
argues that this correction amounts to up to 1.4 mags in the 
B band. Fluctuations in the dust corrections can be of the 
same order of the offset we see between model and data. The 
better agreement between model and data in the K band is 
therefore particularly significant in this scenario, since this 



band is basically unaffected by dust reddening (and in ad- 
dition the fc-correction performed on the data to obtained 
rest-frame magnitudes is relatively small). However, if in- 
ternal extinction corrections were the main cause of the dis- 
agreement in the optical TF, this would imply that in all 
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the observed relations the dust reddening has been sistem- 
atically underestimated, which seems implausible. 

The disagreement between model and data is present in 
all bands, but it is much worse in the optical/blue bands, and 
presents a morphology dependence. This suggests that the 
cause may be found in the balance between the emissions 
of stellar populations of different ages, which means that 
the star formation histories in the model galaxies are not 
realistic. An excess of luminosity like the one we obtain, 
with a variance from B to K, may hint to an excessively 
slow decline of the star formation rate towards low redshifts. 
As will be shown later on, another hint in favour of this 
argument is the redshift evolution of the TF, the fact that 
the match between model and data gets better at higher 
redshifts. 



4.2 Morphology dependence of the TuUy-Fisher 
at z=0 

From the previous Section it is evident that the ratio 
Afbuige/Afdisk plays a role in determining the slope of the 
model TF. In particular, as the bulge mass increases with 
respect to the disk, the slope of the model TF decreases 
and the scatter increases. In general, model late-types agree 
better with the observed slope. This is encouraging, as the 
samples of galaxies selected to observe the TF relation are 
usually biased towards late-type spirals, which tend to fea- 
ture cleaner rotation curves, due to stronger emission lines 
(Masters et al., 2008) 

The Mbuige/Mdisk ratio is a reasonable proxy for Hubble 
type, in that it gives a fair representation of the percentage 
of galaxy mass that is involved in star formation activity 
(model bulges do not receive cooling of fresh gas). Fig. ((5]) 
shows the morphology dependence of the 2 = model TuUy- 
Fisher relation in the K band, according to our bulge/disk 
stellar mass classification. The results are shown in compar- 
ison with data from Masters et al. (2008; dashed black lines 
and yellow shaded area representing the relation and its scat- 
ter), and are divided according to Hubble type. In the up- 
per panel, we compare observed Sa galaxies with the model 
galaxies with 0.5 < Mbuige/Afdisk < 0.75 {red dots), and 
we also show model galaxies with 0.75 < Mbuige/Mdisk < 1 
{red circles); in the middle panel, we compare Sb types with 
model galaxies with 0.25 < Mbuigo/A/disk < 0.5 {red dots); 
in the lower panel we compare Sc types with model galaxies 
with 0.25 < Mbuigc/Afdisk {red dots). 

The model TF shows a differentiation with galaxy mor- 
phology, contributed by the more sophisticated recipe for 
the model rotation curves. The same differentiation is shown 
by the data, with later-type spirals exhibiting a shallower 
slope. The model TF is consistent with the data, with a 
good match with Masters et al. (2008) for Sb/Sc spirals (see 
Table [l]), for velocities log{W) > 2.2, although the slope is 
slightly different and the scatter is smaller. The model re- 
lation for Sa spirals is on the faint side of the data, with a 
clear bending of the slope at the high-mass end, due to the 
increasing bulge/disk mass ratio in the model galaxies. The 
success of the model in producing a morphology dependence 
of the TF shows that the balance between bulge and disk 
emission in the K band is correctly accounted for thanks 
to the use of the M05 stellar populations. Moreover, it also 



shows that the balance between the bulge and disk dynamics 
is effectively represented by the new rotation curve recipe. 

As evident from Fig. ([5]) and the Figures in the previ- 
ous Section, the model TF relation shows a small scatter 
for late-type galaxies, and an increasingly large scatter eis 
the A/buigo/Afdisk ratio increases. This mirrors the different 
mechanism for bulge formation in late and early type spi- 
rals. Small bulges are likely to be formed entirely through 
secular evolution, they do not alter the rotation curve in a 
significant way, and their stellar populations are coeval with 
the ones in the disk. Massive bulges are more likely to be 
formed via mergers, which introduces a large random factor 
both in the ratio Mbuigo/Mdisk and in the relative ages and 
chemical composition between disk and bulge stellar popu- 
lations. The presence of a large bulge manifests itself with 
a steepening or a bump in the central part of the rotation 
curve, which increases the rotational velocity at 2.2 Ru at a 
given magnitude. Moreover, the presence of a massive bulge 
tends to dim the galaxy emission per unit mass (with respect 
to a disk-dominated object of the same mass). In addition 
to increasingly massive bulges, the downward trend at the 
very high-mass end of the theoretical relation, for the earlier 
types, is due to the onset of AGN feedback, that quenches 
star formation, and makes galaxies less luminous per unit 
mass. 

Although we plotted the Masters et al. (2008) relation 
as a straight line, from their paper it is evident that this 
trend is visible in the data as well, and is apparent for Sa 
and Sb types. The model seems to mirror this behaviour, 
although in the data, it is observed at about log{W) ~ 2.6, 
while GallCS produces it around log{W) ~ 2.4 — 2.5. It 
is evident that GallCS cannot cover the mass range of the 
data due in part to the simulation limits; very large spirals 
are extremely rare in nature, and they occur with prob- 
ability in the simulation. Moreover, for increasing masses 
the model galaxies are more likely to experience mergers, 
which produce larger bulges and increase the Mbuigc/Afdisk 
ratio; galaxies therefore are transformed into earlier types 
and progressively disappear from the TF relation. 



4.3 Relative contribution of dynamics and stellar 
population models 

Fig. ^ summarizes the relative effects of a more precise 
and realistic model rotation curve, and of more complete 
stellar populations model. In the four panels, the z — 
model TF relation in the K band is shown for the Sa spirals 
(with errorbars representing its scatter), compared to the 
corresponding observed TF by Masters et al. (2008) {black 
solid line, see lower panel of Fig. ((SI). In each panel, a dif- 
ferent combination of dynamical model — stellar population 
model is shown, between the 'new' and 'spherical' model for 
the rotation curve (see Section 2.1) and between M05 and 
PEGASE stellar population models (see Section 2.2). The ef- 
fect of different rotation curve models and stellar population 
models on the TF is of comparable magnitude. The M05 run 
of the semi-analytic model produces a K-band TF brighter 
by ~ 0.4 mags than the PEGASE run. The 'new' recipe for 
the rotation velocity shifts the TF by ~ 20 — 30% in velocity, 
depending on the galaxy mass and morphology. There is a 
certain degeneracy between the dynamics and stellar pop- 
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Figure 5. The morphology dependence of the 2 = K-band TuUy-Fisher relation. The model predictions are compared with with data 
from Masters et al. (2008) {dashed solid lines, yellow shaded area representing the data scatter). From upper to lower: observed Sa 
galaxies compared with model galaxies with 0.5 < M]-mige/M(jigi5 < 0.75 {red dots) and galaxies with 0.75 < Mbuigo/Mdiaj; < 1 {red 
circles); observed Sb galaxies compared with model galaxies with 0.25 < Mbuigc/Af^isk < 0.5 {red dots); observed Sc galaxies compared 
with models galaxies with 0.25 < Mbuigp/M^isk ('"^'^ dots). 



ulation models. The 'new' dynamical model works better 
than the 'spherical' for both stellar population models. 

The analysis carried out so far implies that the model 
TF is flawed at z = 0, in the B, I and K bands, even after 
we adopted refined and updated prescriptions for the galaxy 
rotation curves and the stellar populations. We use this con- 
clusion to gain a better grasp of more profound problems of 
the model in reproducing the spiral population, in particu- 
lar regarding the cooling and star formation histories at low 
redshifts, as will be discussed further in this work. 

4.4 Satellites and the TF relation 

Satellite galaxies in the simulation are defined as objects 
that do not reside in the center of the parent dark mat- 
ter halo, but are orbiting around or infalling into a central 
galaxy. The gas in satellites is stripped because of ram pres- 
sure and tidal effects, while the cooling of new gas is pre- 
vented. This effectively shuts down star formation in satellite 
gala:xies, which evolve passively unless a merger event takes 
place. For this reason, model satellite galaxies are ill-suited 



to be compared with real spirals. For the purpose of clarifi- 
cation and completeness, we show the predicted K-band TF 
relation for satellite spirals and compare them to centrals. 

Fig. (O shows the z = K-band TF relation for a sam- 
ple of model spiral galaxies in the range < Mbuige/Mdisk < 
1, and additionally split into central {magenta dots) and 
satellite {blue dots) galaxies. It is clear that satellites do not 
follow a TF relation, but scatter below the central galaxy 
relation. The main reason is that satellites are less bright, 
due to the quenching of the star formation which, even if 
not affecting directly the K-band, raises the overall mass- 
to-light ratio. This affects expecially the low mass end of 
the satellite distribution. 

Note that a significant fraction of the galaxies in the 
Master et al. samples are cluster galaxies, and they follow 
the same TF relation as the field samples. We can there- 
fore use the TF relation as a tool to evaluate the ability 
of the model to reproduce the structure of satellite galax- 
ies. We conclude without doubt that the model cannot pro- 
duce realistic satellites. Again, the main problem resides in 
the star formation history. Real satellite spiral galaxies, al- 
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Figure 6. Relative contributions of the rotation curve recipe described in Section 2.1 and the stellar population models described 
in Section 2.2. In the four panels, the z = model TF relation in the K band is shown for the model Sa spirals, compared to the 
corresponding observed TF by Masters et al. (2008) {black solid line, see lower panel of Fig. JSJ). In each panel, a different combination 
of dynamical model — stellar population model is shown, between the 'new' and 'spherical' model for the rotation curve (see Section 2.1) 
and between M05 and PEGASE stellar population models (see Section 2.2), as indicated by the labels in each panel. 



though perturbed when in dense environments, are still star- 
forming, in particular those selected for TF studies. Hierar- 
chical models, in order to reproduce the observed colour- 
magnitude relation, are forced to shut down star formation 
in satellites. This highlights a fundamental problem in hi- 
erarchical models, that certainly deserves further investiga- 
tion. 



5 STELLAR MASS VS K-BAND TF 
RELATION: SO GALAXIES? 

There is an ongoing debate about the formation of SO galax- 
ies. The origin of these objects is still not understood, 
but current scenarios include massive bulge formation, tidal 
stripping and star formation quenching as possible causes of 
a transition from spirals to SO, via a progressive fading of 
the disk (see Poggianti et al. 2001 and references therein). 
In particular, the possibility that SO galaxies may be spi- 
ral galaxies where star formation has been shut down, and 
the spiral structure has thus become invisible, can be tested 



with the present study. An interesting consequence of this 
is the fact that the TF relation should present an offset be- 
tween SO and spiral galaxies (Williams et al. 2009, Bedregal 
et al. 2006, and references therein), which should be visible 
in the photometric TF but not in the stellar mass-velocity 
TF relation. However, Williams et al. (2009) find an offset 
between observed spirals and SO both in the K-band TF and 
in the stellar TF, and therefore argue that SO galaxies are 
not fading spirals. We use our model to verify whether star 
formation quenching introduces an offset between the stellar 
and K-band TF relation. 

In this Section we show the model TF relation in the 
K-band, compared to the stellar TF relation Mstar vs Kot, 
for star-forming and nearly-passive disk galaxies at z = 
0. We again differentiate the galaxy morphology through 
the Mbuige/Mdisk ratio. Model disk galaxies with large 
Afbuigo/A^disk and low star formation rate should in prin- 
ciple present a spectral energy distribution similar to real 
SO galaxies. First, the large bulges (built-up preferentially 
via minor mergers) and the small gaseous mass make these 
objects the transition point between spirals and ellipticals 
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Figure 7. The predicted z = K-band TuUy-Fisher relation, for a sample of model spirals in the range < Mbuige/Mdisk < Ij and 
additionally split into central {magenta dots) and satellite (blue dots) galaxies. The model is compared with the observed relation by 
Masters et al. (2008; solid black line for the relation, and dot-dashed lines for its scatter) and Verheijen 2001 {green solid/dashed lines). 



in the model. Second, the low star formation rates mimic 
the quenching due to gas stripping in dense environments 
as well as the natural fading due to the exhaustion of the 
gas reservoir. 

Fig. ((8| shows the comparison between the z — stellar 
TF {left panels) and K-band TF {right panels) relations, for 
two subsamples of model galaxies divided according to the 
instantaneous SFR. In the upper panels, galaxies are selected 
depending on their total star-formation rate, and triangles 
show star-forming galaxies with SFR > 3 Mq/jt, while cir- 
cles represent nearly-passive galaxies with SFR < 1 Mq /yr. 
In the lower panels instead, galaxies are selected according 
to their specific star formation rates, and triangles represent 
star-forming galaxies with sSFR > 10"^" yr~^, while circles 
represent nearly-passive galaxies with sSFR < 3 10~^^ yr~^ 
(the separation between star-forming and nearly-passive was 
taken from Milky- Way values as a reference, see Munoz- 
Mateos et al., 2007). In all panels, galaxies are also colour- 
coded according to morphology, with cyan circles/filled tri- 
angles representing 'late-type' spirals with Mbuige/Afdisk < 
0.5, and red circles/filled triangles representing 'early-type' 
spirals with 0.5 < A/buigc/Afdisk < 1- Red circles therefore 
represent the model rendition of SO-like objects (in terms of 
stellar populations). 

When the galaxies are selected according to the total 
star-formation rate, we find a continuous stellar TF relation, 
as expected, between star-forming and nearly-passive galax- 
ies, while we find a double sequence in the K-band, with the 
star-forming galaxies significantly brighter than their pas- 
sive counterparts, for a given mass (the difference in mag is 
~ 0.5, consistent with the findings of Williams et al.). We 
also find a morphology segregation in the K-band TF, with 



later-type spirals at the brigher end of the population (ex- 
cept the very top of the mass distribution, where later-types 
are not found in the model). This segregation tends to dis- 
appear in the stellar TF. This result seems to confirm that, 
indeed, SO can be considered quiescent disks that are in the 
process of shutting down the star formation, and in this case 
the model SO are represented by the red circles. 

If we select the galaxies according to the specific star- 
formation rate, we find again a different trend between the 
stellar and K-band TF relations. The stellar TF relation 
shows a net segregation between star-forming and nearly- 
passive galaxies, which form a separate sequence. In addi- 
tion, there is a net morphology segregation, which puts the 
earlier- types, nearly-passive galaxies at the massive end of 
the relation for any given velocity. In the K-band relation the 
star-forming, later-type galaxies are clearly on the brighter, 
lower-mass end of the relation for any given velocity. If we 
consider that the red circles as the model SO galaxies, the 
picture is consistent with SO galaxies being fading spirals. 
Earlier types, characterized by massive bulge formation, oc- 
cupy the high-mass end of the population. The presence of 
a massive bulge reduces the fraction of galaxy mass that 
is actively star-forming, so that the specific star-formation 
rate and the luminosity are low compared with later-types 
of similar mass. 

Although more investigation is needed on this topic, we 
are encouraged to conclude that the scenario that describes 
SO galaxies are fading spirals is promising. We note that a 
direct comparison between our results and the findings of 
Williams et al. (2009) is problematic, due to the differences 
between the modeling of the dynamical mass by Williams et 
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Figure 8. Comparison between the stellar and K-band TF relations for star-forming and non-star forming model galaxies. Left panels: 
model stellar TF relation at z = 0. Right panels: model K-band TF relation at 2 = 0. Upper panels: galaxies selected depending on their 
total star- formation rate. Triangles represent star forming galaxies with SFR > 3 MQ/yr and circles represent nearly-passive galaxies 
with SFR < 1 M0/yr. Lower panels: galaxies selected depending on ther specific star-formation rate. Triangles represent star forming 
galaxies with sSFR > 10"^" yr~^ and circles represent nearly-passive galaxies with sSFR < 3 10~^^ yr~^- In all panels, galaxies 
are colour-coded according to morphology, with cyan circles/filled triangles representing 'later-types' with Mbuige/Mdisk < 0-5: a^nd '"^ci 
circles/filled triangles representing 'earlier-types' with 0.5 < M^^^i^^/M^^^]^ < 1. 



al. and the mass- velocity relation in GallCS. Further inves- 
tigations on these results are under development 



6.1 Evolution of the optical TF 



6 THE REDSHIFT EVOLUTION OF THE TF 
UP TO Z=l 

The evolution of the TF is a good test of the galaxy assembly 
mechanism in the semi-analytic model, and here we present 
it in the optical and near-IR. 



Fig. ((9)1 shows the redshift evolution of the model rest-frame 
B-band TF relation from z ~ 0.4 to 2 ~ 0.8, compared with 
data from Fernandez-Lorenzo et al. (2009; black solid/dashed 
lines representing the observed relation and its scatter, in all 
panels) , who derive the galaxy rotational velocity in a some- 
what different manner than the rest of the published obser- 
vational relations, based on optical line widths (see also Ver- 
heijen 2001 for a comparison between TF relations obtained 
with different velocity estimators). In the z — 0.4 panel, data 
from Bohm & Ziegler (2007) are shown as black triangles, 
and the relation by Bamford et al. (2006) at < 2; >= 0.33 is 
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Figure 9. The redshift evolution of the rest-frame B-band TuUy-Fisher relation from z ~ 0.4 to z ~ 0.8, eompared with data from 
Fernandez-Lorenzo et al. (2009). Solid/dashed black lines represent the observed relation and its 2a uncertainty, in all panels. In the 
z = 0.4 panel, the data from Bohm & Ziegler (2007) are shown as black triangles. The relation by Bamford et al. (2006) at < z >= 0.33 
is also shown (red line). The model galaxies are represented by dots, colour-coded as in the previous Figures. 



also shown for comparison. The colour-coding of the model 
dots is as in the previous Figures. 

This figure shows a very good agreement between the 
model late-type spirals (Mbuigc/Afdisk < 0.5, yellow and ma- 
genta points) and the data at 2 ~ 0.6 and z ~ 0.8, in slope, 
zero point and scatter. At z ~ 0.4 the model reproduces the 
observed slope, but tends to be slighly fainter in zero point 
(although still inside the observed 2a uncertainty). Again, 
the model galaxies show a very clear morphological differ- 
entiation, and the earlier-types fall off the relation at the 
high-mass end at all redshifts, which as already discussed, is 
an effect mainly due to bulge formation, and AGN feedback. 

The decline of the global SFR with time, together with 
the migration of the star formation to lower and lower galaxy 
masses (the downsizing effect), causes the galaxies to be- 
come fainter and the slope of the TP to get shallower at 
lower redshifts (also confirmed by the results of Weiner et 
al. 2006). 

The model optical TF agrees very well with the data at 
high redshift, and by z ~ 0.4 the match starts to fail, where 
the model galaxies are too faint. By z ~ 0, the model is 
instead too bright (Fig. Q. This seems to confirm our sus- 
picion that the model galaxy star formation histories are the 
main drivers of the discrepancy (or match) between model 
and data. At z > 0.5 the model spiral gala^xies seem to 
be quite realistic, while for lower redshifts the gas cooling, 
the accretion of substructures and the supernovae feedback 
conspire to produce an excessive star formation at odds with 
observations. 

6.2 Evolution of the near-IR TF relation 

Fig. (jlOp shows the evolution of the observed-frame K-band 
TF relation, in the redshift range from z ~ 0.2 to z ~ 1.2. 
The colour-coding of the model galaxies is the same as in the 
previous Figures. We compare the model predictions with 
data from Fernandez- Lorenzo et al. (2010; solid black trian- 
gles) and Bohm & Ziegler (2007; solid black squares). In this 
case the observed magnitudes are not corrected for dust ex- 



tinction, so the model galaxy spectra are reddened according 
to their star-formation rate, as discussed in Section 2. 

The scatter in the data here is large, and we do not 
attempt to determine the observed slopes and zero points. 
Nonetheless, the accord between the model TF and the ob- 
servations is good in the range 0.4 ^ z ^ 1.2. Both the 
observed velocity range and luminosity range are well re- 
produced. The model galaxies also occupy the same velocity- 
luminosity space as the data. Moreover, the redshift evolu- 
tion is reasonably well reproduced. 

At z ~ 0.2, the model galaxies are too bright at the 
low- mass end, consistently with the results at z = (Fig.|4]). 
Following the combined information of the evolution of the 
TF in the optical and near-IR, we argue that the model star 
formation histories are realistic down to z > 0.4, and odd 
for lower redshifts. 



7 DISCUSSION 

Disentangling the contributions of the different physical pro- 
cesses conspiring to produce the TuUy-Fisher is not a trivial 
task. The slope and zero-point of the relation are determined 
by the interplay of the dynamics of galaxy assembly inside 
the dark matter halo, the cooling that regulates the gaseous 
and stellar content in each halo, the mass accretion and star 
formation histories and the supernovae feedback, all of which 
affect the structure and stability of disks as well as the evo- 
lution of the luminosity across all bands. Thus, matching the 
observed TF scaling relation is a fundamental check that the 
model is correctly predicting the formation and evolution of 
spirals. 

However, a meaningful comparison with the observed 
TF necessitates of some fundamental steps, that ensure that 
there are no systematics that can offset the results. The first 
step is the correct modeling of the galaxy rotation curves, 
which are not directly predicted by the semi-analytic model 
itself, but need to be determined based on the galaxy fun- 
damental dynamical quantities. The rotation curves extrap- 
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Figure 10. The redshift evolution of the observed-frame K-band TuUy-Fisher relation from z ~ 0.2 to z ~ 1.2, compared with data from 
Fernandez-Lorenzo et al. (2010; solid black triangles) and Bohm & Ziegler (2007; solid black squares). The model galaxies are represented 
by dots, colour-coded as in the previous Figures. 



olate spatial information in model galaxies, following a dy- 
namical recipe that is physically motivated and that needs 
to be as sophisticated as possible. In our 'new' model, we 
made use of models for the mass distribution of the galactic 
components and angular momentum transfer to the disk ac- 
companying bulge formation. It is worth noting that a com- 
plementary tool in this sense is represented by SPH simula- 
tions, which favour a detailed study of the galaxy structure, 
and are able to reproduce the TF relation for single galax- 
ies, (e.g. Piontek & Steinmetz, 2009b), even if the shape 
of the rotation curve often remains unrealistic (see for in- 
stance Governato et al. 2007), and massive spiral galaxies 
still represent a problem. The sophistication of current SPH 
simulations also allows for detailed studies of the effects of 
perturbations such as ram-pressure on the rotation curves 
and the star formation (Kapferer et al. 2009, Kronenberg 
et al. 2008). However, the current resolution limits for this 
kind of simulations does not allow for cosmological runs and 
a statistical study. 

The second step is the implementation of comprehensive 
stellar population models, which assures that the photomet- 
ric side of the scaling relation is not marred by systematic 
offsets. We showed that this element can introduce a bias 
as large as 1 mag in the K-band, and 0.3-0.4 mag in the I 
band. 

The third step is a consistency check in the compari- 



son with observations, that presents challenges not related 
to the semi-analytic model itself. For instance, in a given 
band, observed TF relations sometimes are not consistent 
with one another, and the reason for these discrepancies is 
probably to be found in the data analysis. In fact, a sig- 
nificant degree of modeling of the raw data is required to 
produce an observed TF. To name a few, corrections are 
required for completeness, internal extinction, inclination, 
average distance and size in case of clusters, morphology, 
peculiar velocities (see for instance Masters et al. 2006). In 
addition, observed-frame magnitudes have to be converted 
into rest-frame through fc-correction. and this introduces a 
further bias in the comparison. As we pointed out in Tonini 
et al. (2010), it is preferable to translate the semi-analytic 
model into observed-frame and compare the predictions with 
the apparent magnitudes. This is crucial expecially for star- 
forming objects, because their evolution is fast and the cor- 
rections are significant, (expecially for internal extinction 
and fc-correction) . 



After the systematics in the modeling and in the com- 
parison with data are taken care of in 'post-production', the 
hierarchical model can be truly tested, and the problems of 
the model in reproducing the TF relation can be investi- 
gated. 
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T.l Star- formation histories 

The offset between the model and the observed TF relation 
has the following features: i) it depends on redshift, with 
a better agreement between model and data at higher red- 
shift; a) at redshift 2 = it depends on the photometric 
band, with a better agreement in the K band and the worst 
performance of the model in the B band; in) it presents a 
morphology dependence across all redshifts and photometric 
bands. These facts suggest that the star formation history of 
the model galaxies is unrealistic, after z ~ 0.4. As a conse- 
quence, the predicted aging of the stellar populations is off, 
and the emission of stellar populations of different ages is un- 
balanced. In particular, the decline of star formation with 
time seems not to be fast enough, with the consequence that 
the fraction of the total stellar mass constituted of young 
and intermediate-age stellar populations is too high, a fact 
that predominantly affects the UV and optical bands, and 
the near-IR via the TP-AGB light. Note that this does not 
imply a wrong instantaneous SFR at z = 0, which in fact is 
predicted to be consistent with observations (Hatton et al. 
2003). 

Central galaxies in a hierarchical model never stop ac- 
creting material, including fresh gas and gas previously ex- 
pelled from the galaxy itself. Although this type of star for- 
mation history can be chaotic, in most cases at low red- 
shifts GallCS star formation histories show a sufficiently 
steady decline to be approximated in the form of r— models: 
SFR{t) = SFRoexp{—t/to). We calculated the predicted 
luminosity of model galaxies with different star formation 
timescales, and found the latter to be a fundamental pa- 
rameter regulating the relative luminosities in B and K. For 
instance, a galaxy of age of 10 Gyr with a star-formation 
history like SFR{t) = SFRoexp{-t/ 10) (with a e-folding 
time r = 10 Gyr) produces a B-band magnitude brighter 
by ~ 0.83 mags compared to a model like SFR{t) = 
SFRoexp{—t/3) (with a e-folding time r = 3 Gyr), while 
it is brighter in the K-band by ~ 0.46 mags. From Fig. 3] it 
is evident that a shift of ~ 0.83 mags in the B band would 
produce a much better agreement between model and data, 
while at the same time, a shift of ~ 0.46 mags in the K band 
would still preserve the accord with the observed relations. 

A faster decline of the star formation in time allows the 
galaxies to fade more rapidly at low redshifts, and brings the 
predicted TF relation down in luminosity. To size this effect, 
consider what happens to satellites after the star formation 
is shut down (Fig.[7|. Also, as shown in all previous Figures, 
as more and more of the fraction of galaxy mass is in the 
inactive bulge component, or AGN feedback kicks in, the 
galaxy tends to fade and fall below the TF relation. 

We argue that, in order to fix this problem, a more rad- 
ical revision of gas infall and cooling is probably needed. A 
more detailed investigation of this point is of great interest, 
but it is beyond the present scope, in that it requires an 
extended comparison with SFR derived from observations 
at different redshifts. Note that the derivation of SFR from 
observations is affected by the adopted stellar population 
models, dust extinction corrections, star-formation history 
laws, as discussed in Maraston et al. (2010), and the use of 
inhomogeneous sets of such priors in comparison with our 
model will lead to significant biases. This line of investiga- 
tion is currently under development 
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Figure 11. Model disk scale-lenghts as a function of I-band lu- 
minosity, compared with the relation obtained by Courteau et 
al. (2007) from observations {solid/dashed lines for the relation 
and its scatter). The empty dots represent model galaxies in the 
'spheroidal' model, colour-coded according to Mbuigc/Af^isk a-s in 
the previous Figures. Filled triangles represent the same galaxies 
after angular momentum transfer from the bulge is evaluated (the 
'new' model). 

In addition to the star formation history, a second-order 
effect on the optical TF is caused by the chemical evolution 
model in use. The current version of GallCS implements the 
Pipino et al. (2009) recipe which, compared to Hatton et 
al. (2003), produces slightly lower total metallicities, caus- 
ing the galaxies to be brighter. Depending on age and star 
formation history, the difference in magnitude amounts to 
anything between and ~ 0.5 in the B band, and less for the 
I and K bands. Accounting for uncertainties in the chemi- 
cal evolution recipe would only lead to a modest increase of 
the TF scatter, which would be consistent with the scatter 
in the data, but it would not shift the overall relation in a 
significant way. 

7.2 Disk scale-lenghts and the disk-halo 
connection 

There is also the possibility that part of the discrepancies 
that we encountered in our comparison with data are due 
to a flaw in the modeling of the galactic dynamics. In par- 
ticular, the galaxy rotation may be offset, regardless of the 
galaxy type, due to a wrong determination of the galaxy an- 
gular momentum. This can be investigate by analyzing the 
disk scale-lenghts Ro produced by the model. 

Figure [11] shows the model disk scale-lenghts as a 
function of I-band luminosity, compared with the rela- 
tion obtained by Courteau et al. (2007) from observations 
(solid/dashed lines for the relation and its scatter). The 
empty dots represent model galaxies in the 'spheroidal' 
model, colour-coded according to Mbuigo/Mdisk as in the pre- 
vious Figures. Filled triangles represent the same galaxies 
after angular momentum transfer from the bulge is evalu- 
ated (the 'new' model). 

It is evident that the model cannot reproduce the 
Courteau et al. relation. The disk-scale lenght appears to 
be too large, and the scatter too big, the latter a problem 
that was mentioned already in Hatton et al. (2003). More 
in detail, the model does not reproduce the Courteau et 
al. relation for the late-type spirals {yellow points), which 
are mostly untouched by our 'new' velocity model. In the 
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old 'spheroidal' model, the earlier-type spirals present a 
marginal agreement with the data, which we argue is ac- 
cidental. In fact, the values of the disk-scale lenghts inferred 
from observations are commonly obtained by fitting the sur- 
face brightness profile with an exponential law, and that the 
result is more accurate the closer the galaxy is to a pure disk. 
Galaxies with large bulges are not ideal candidates, and they 
tend to be excluded from the observational samples. The 
'new' model effectively corrects for the presence of the bulge, 
moving the earlier-type spirals in the same -Rd — Li space as 
the later-types, in a sense making all galaxies consistently 
offset with the Courteau et al. relation. 

Before analyzing the possible causes of the discrepancy, 
it should be noted that a source of error stems from the use of 
Eq. ((3]) for the determination of Rn- This formula, adopted 
by Hatton et al. 2003, produces the highest Rd under angu- 
lar momentum conservation, in the MMW formalism. Ad- 
ditional shape factors that can reduce the value of Rd are 
set to unity, but as discussed in Tonini et al. (2006a), their 
values as inferred from observed rotation curves are actually 
close to unity, and vary slowly with disk mass. 

As for the cause of the discrepancy with the Courteau 
data, we offer 2 explanations: 

I) as discussed in Tonini et al. (2006a) and D'Onghia & 
Burkert (2004), the halo spin parameter, which is found to 
vary from < A >= 0.04 in GallCS to < A >= 0.05 in other 
simulations, is too high compared with the one that can be 
inferred from observations of spirals. If < A 0.03 for in- 
stance (the value inferred by Tonini et al. 2006a) , then Rd is 
smaller by ~ 25%, which would push our model _Rd towards 
the Courteau et al. relation. Note that, although A originates 
naturally from tidal forces during structure formation, the 
discrepancy between the simulated and observed A is not 
necessarily due to a flaw in the simulations, but rather it 
hints at a subsequent redistribution of angular momentum 
between baryons and dark matter; 

II) as mentioned in Section (2.1), the ratio ju/mo is 
equal to unity only under angular momentum conservation, 
which is not a realistic scenario during galaxy formation. 
In particular, if the baryonic component that collapses to 
form the protogalaxy is clumpy, it likely dissipates angu- 
lar momentum in its path to the centre of the dark matter 
halo, as proposed by Tonini et al. (2006b), and therefore 
ju/mu < 1. Consequently, the disk scale-length is smaller, 
and the tension between model Rd and observations is par- 
tially alleviated (jd/'tid ~ 0.75 leads to a decrease of Ru 
by ~ 25%). The MMW recipe can therefore be considered 
a 'maximal scale-length scenario', with the real scale-length 
distribution being centered on smaller values of Rd because 
of angular momentum transfer. 

Both the effects described above contribute in reducing 
the size of Rd, and are currently not investigated in semi- 
analytic models. Note also that such effects do not depend 
on redshift nor on photometric band, therefore they do not 
affect our conclusions on the model star formation histories 
discussed in the previous section. 

An additional, largely unknown source of scatter, on 
both the disk scale-lengths and the rotation curve, is rep- 
resented by the dark matter halo structural evolution. On 
the one side, readjustments in the halo profile induced by 
the baryon infall can lead to a dramatic change in the in- 
ner density profile (Tonini et al. 2006b), thus altering the 



contribution of the halo to the total velocity. On the other 
side, tidal interactions between halos and mass accretion in 
the hierarchical assembly repeatedly strip and add material 
in each halo. There is no standard recipe in semi-analytic 
models to recalculate the halo equilibrium structure that 
takes into account all these effects, and since there is no 
spatial information inside single objects, the mass is the pri- 
mary driver of halo evolution. The virial radius in particular 
is not recalculated after every interaction, with the conse- 
quence that the halo density profile can oscillate quite a lot. 
Note that also in N-body simulations, although all halos are 
fitted to a common functional profile (like the NFW or the 
isothermal sphere), the scatter on the profile is very large, as 
is the scatter on the mass-virial radius relation (Bullock et 
al. 2001). The scatter in the virial radius directly translates 
into a scatter in Rd, while the scatter in the density profile 
affects the halo velocity. In particular, an overestimation of 
the disk scale-lengths will occurr in cases where the mass 
loss via tidal stripping dominates the halo evolution. 

The offset between model and observed disk scale- 
lengths represents an interesting failure of the model, that 
arguably sheds light on some missing physics. These con- 
siderations suggest that baryon cooling and halo structural 
evolution in the model need some more fundamental revi- 
sion, which is beyond the scope of the present work. 



8 SUMMARY AND CONCLUSIONS 

We analyzed the predictions of the hierarchical galaxy for- 
mation model GallCS on the TuUy-Fisher relation and its 
evolution with redshift. We introduced two new elements 
in the model: 1) a new recipe for the determination of the 
galaxy rotation curves, based on the correct dynamical mod- 
eling of the different galactic components, and taking into 
account the angular momentum exchange between bulge and 
disk due to secular evolution; 2) the M05 stellar popula- 
tion models, that include an exhaustive treatment of the 
TP-AGB emission, which affects galaxies with ongoing star 
formation. 

Our main conclusions are: 

• the improvement on the dynamical description of disk 
galaxies and on the stellar population models impacts the 
model TF by effects of comparable magnitude; in particular, 
velocities are shifted to values lower by ~ — 50% depend- 
ing on the galaxy mass and morphology, and the K-band 
magnitudes are brighter by ~ 0.2 mags at redshift zero and 
up to > 1 mags at redshift z ^ 1 and above; 

• at redshift z — the model reproduces the slope of the 
observed TF for Sb/Sc spirals, in the B, I and K band, but 
not the zero-point, which is too bright in all bands. In par- 
ticular, in the K band the zero-point is too bright for Sb/Sc 
galaxies at the low-mass end {log{W) < 2.3), although the 
model galaxies lie within the scatter of the observed TF rela- 
tion. In the B and I band, the predicted zero-point for Sb/Sc 
galaxies is too bright across the mass range. We argue that 
the most probable cause of the discrepancy lies in unrealistic 
star-formation histories at z < 0.4; 

• the model predicts a morphology dependence of the 
TF relation at all redshifts. At z = observations in the K- 
band confirm this trend, and the comparison between model 
and data is encouragingly good, expecially for Sb/Sc spirals. 
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This is an important confirmation that our more sophis- 
ticated treatment of the galaxy dynamics is reahstic. The 
model rotation curve now mirrors the morphological differ- 
entiation observed in nature, well representing the balance 
between bulge and disk, and the consequent variations in the 
predicted TF along the Hubble sequence are put in sharper 
evidence. The model predicts a steeper slope and a smaller 
scatter for later-type galaxies, and a progressive flattening 
of the slope and larger scatter for earlier types. This is co- 
herent with the picture of bulge formation in a hierarchical 
scenario, where the incidence of mergers is significant in the 
case of massive bulges, while small bulges are more likely to 
develop via secular evolution. Moreover, with the M05 stel- 
lar population models, GallCS is able to correctly reproduce 
the balance between the bulge and disk emission in the K 
band; 

• the scatter in the model galaxies increases from late 

types (Sb/Sc) to early types (Sa/SO), because the main 
mechanism for bulge formation switches from secular evo- 
lution (Sb/Sc) to minor mergers (Sa/SO). At the high-mass 
end, earlier-type galaxies show a flattening of the slope and 
tend to fall below the relation due to the presence of massive 
bulges and AGN feedback; 

• the study of the TF relation for model satellites galax- 
ies confirms previous hints that the model cannot produce 
realistic non-central spirals; 

• the TF relation can be used as a discriminating tool 
to investigate the origin of SO galaxies, in the scenario where 
these objects form from the fading of spirals after star- 
formation shutdown; 

• the model reproduces the rcdshift evolution of the 
TF for later-type galaxies from z = 0.8 to z = 0.4 in the 
rest-frame B band, where slope, zero-point and scatter are 
reproduced for Sb/Sc spirals. The model is consistent with 
the data in the observed-frame K band at redshifts between 
z = 1.2 and z — 0.4. The simultaneous match in the optical 
and near-lR evolution, although needing more tests with 
larger data samples, indicates that the mass assembly and 
star formation histories, as well as the supernovae feedback 
implementation, are satisfactory at these redshifts; 

• the model cannot reproduce the whole optical Tully- 
Fisher evolution from 2~0toz~1.2. We argue that the 
decline of the star formation is too slow, thus preventing the 
model from matching simultaneously the local and distant 
relations; 

• the model produces disk scale-lengths too large com- 
pared with observations. We argue that this mirrors a too 
simplistic recipe of disk formation inside dark matter halos, 
that does not take into account angular momentum redis- 
tribution during the baryonic collapse. 
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